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ABSTRACT 

We study here the formation of heavy r-process nuclei in the high-entropy environment of rapidly expanding neutrino- 
driven winds from compact objects. In particular, we explore the sensitivity of the element creation in the A J; 130 
region to the low-temperature behavior of the outflows. For this purpose we employ a simplified model of the dynamics 
and of the thermodynamical evolution for radiation dominated, adiabatic outflows. It consists of a first stage of fast, 
exponential cooling with timescale Td yn , followed by a second phase of slower evolution, assuming either constant density 
and temperature or a power-law decay of these quantities. These cases describe a strong deceleration or decreasing 
acceleration of the transsonic outflows, respectively, and thus are supposed to capture the most relevant effects associated 
with a change of the wind expansion behavior at large radii, for example because of the collision with the slower, 
preceding supernova ejecta and the possible presence of a wind termination shock. We find that for given entropy, 
expansion timescale, and proton-to-baryon ratio not only the transition temperature between the two expansion phases 
can make a big difference in the formation of the platinum peak, but also the detailed cooling law during the later phase. 
Because the nuclear photodisintegration rates between about 2 x 10 s K and roughly 10 9 K are more sensitive to the 
temperature than the neutron-capture rates are to the free neutron density, a faster cooling but continuing high neutron 
density in this temperature regime allow the r-process path to move closer to the neutron-drip line. With low (7,n)- 
but high /3-decay rates, the r-processing does then not proceed through a (7,n)-(n,7) equilibrium but through a quasi- 
equilibrium of (n,7)-reactions and /3-decays, as recently also pointed out by Wanajo. Unless the transition temperature 
and corresponding (free neutron) density become too small (T ^ 2 x 10 8 K), a lower temperature or faster temperature 
decline during the slow, late evolution phase therefore allow for a stronger appearance of the third abundance peak. 

Key words. Nuclear reactions, nucleosynthesis, abundances — Stars: supernovae: general — Stars: winds, outflows — 
Stars: neutron 



1. Introduction 

Heavy nuclei beyond the iron peak are known to be pro- 
duced in nature mainly through neutron capture reactions 
(Burbidge et al. 1957). Rapid neutron capture and the re- 
verse photodisintegration processes achieve an equilibrium 
among the isotopes of each heavy element. Then beta-decay 
occurs that leads to the increase of the nuclear charge and 
formation of a new element. When the neutron capture rate 
is much higher than the beta-decay rate (A n7 > Xp) and 
Tg ~ 1 (Tg is defined as the temperature normalized to 
10 9 K), the r-process can start and, for a sufficiently high 
neutron-to-seed ratio, the wave of nucleosynthesis drives 
the process to heavier nuclei, forming, in part, the abun- 
dance r-process peaks at ^80, 130 and 196. 

Although the r-process sites remain unknown, many as- 
trophysical models and sources for r-process elements have 
been proposed during the last 50 years, including, in par- 
ticular, such scenarios as that of an explosion on a neutron- 
star surface (Bisnovatyi-Kogan and Chechetkin 1979), a 
collision of a neutron star with a black hole (Lattimer & 
Schramm, 1976), an explosion of a low-mass neutron star 
(Imshennik, 1992), and the hypothetical escape of nucleon 
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bubbles in case of a soundless stellar collapse (Imshennik & 
Litvinova 2006). The latter scenario was suggested as origin 
of gamma-ray bursts but could also be an interesting site 
for the nucleosynthesis of heavy elements. 

Presently, however, it seems most likely that rapid 
n-capture nucleosynthesis can take place during differ- 
ent stages of supernova explosions (e.g., Hillebrandt 1978; 
Woosley & Hoffman 1992), or in neutron star mergers (e.g., 
Lattimer & Schramm 1974; Symbalisty & Schramm 1982; 
Freiburghaus et al. 1999). An overview of the currently dis- 
cussed possible sites is given in a recent review paper by 
Arnould et al. (2007). 

Supernovae and neutron star mergers have different ad- 
vantages and weak points, but the main difference lies prob- 
ably in the initial neutron-to-proton ratio, which is neces- 
sary to support a sufficiently high free neutron density dur- 
ing several hundreds of milliseconds. With respect to this 
parameter, all astrophysical scenarios can be separated into 
two distinct groups, in which nucleosynthesis can be carried 
out over short or long timescales (so-called short-time or 
long-time solutions, respectively), as introduced by Seeger 
et al. (1965). In nature, the r-process might be realized 
in both these types of sources and also other ones might 
contribute (see for example Cameron, 2003; Arnould et al. 
2007) . It is only future modelling and observations that will 
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be able clarify the multiplicity of the r-process models and 
types, as well as the contributions from different produc- 
tion sites to the r-process element abundances observed in 
the Solar System. 

Among the astrophysical events proposed as sites for 
the r-process, supernova explosions still remain the prefer- 
able ones (see, e.g., Wanajo & Ishimaru 2006, and ref- 
erences therein). In particular, supernova explosions can 
distribute r-process material all over the Galaxy, and esti- 
mated amounts of heavy elements produced in SN explo- 
sions are in accordance with the observations. 

The neutrino-driven wind from a hot neutron star pro- 
duced in a supernova explosion has been considered as a 
probable site for the r-process by many authors (see e.g. 
Meyer et al. 1992; Woosley et al. 1994; Witti et al. 1993; 
Otsuki et al. 2000; Sumiyoshi et al. 2000; Terasawa et al. 
2001; Wanajo et al. 2001). A part of the surface material 
of a neutron star is heated by the supernova neutrinos and 
gets ejected. It can be described as a hot outflow with a 
fairly high entropy and a moderate density. 

Various studies of the r-process element formation in 
proto-neutron star winds have been conducted during the 
past years. Hoffman, Woosley & Qian (1997), employing 
the analytic wind model of Qian & Woosley (1996), ex- 
plored the possibility of third peak creation for different 
combinations of the determining parameters of Y c (the 
elcctron-to-baryon ratio) , entropy, and expansion timescale. 
Altogether, they showed (see their Fig. 10) that for typical 
values of Y c ^ 0.4 in the wind the third r-process peak can 
be produced for combinations ranging from moderate en- 
tropy (s ~ 100 in units of Boltzmann's constant per baryon) 
and very short timescales (t cxp a few milliseconds) to high 
entropies (s <; 400) and long expansion timescales (a few 
100 milliseconds). The wind models existing at that time 
failed to provide the necessary conditions (Witti et al. 1994, 
Qian & Woosley 1996). 

Subsequently, several other studies of neutrino-driven 
winds in the framework of general relativistic gravity 
(Otsuki et al. 2000, Sumiyoshi et al. 2000, Thompson et 
al. 2001, Wanajo et al. 2001, 2002) confirmed the need of 
fairly extreme conditions concerning expansion timescale 
or entropy for strong r-processing up to A ~ 200. Since 
very high entropies could not be obtained in the wind sce- 
nario or were found to be associated with too low mass 
loss rates for any significant production of r-nuclei, these 
studies also demonstrated a preference for the case of mod- 
erate entropies, s ~ 100-200, and very short timescales. 
This seemed to give a bias for winds from compact neutron 
stars with a large mass, M ^ 2M Q , and a small radius 
R ^ 10 km as the most likely site for r-process element 
formation up to the platinum peak. 

This conclusion was found to hold independent of 
whether freely expanding, transsonic "wind" outflows were 
considered (Hoffman et al. 1997, Thompson et al. 2001) or 
subsonic "breeze" solutions (Otsuki et al. 2000; Sumiyoshi 
et al. 2000; Wanajo et al. 2001; Terasawa et al. 2001, 2002) 
or supersonic winds with some fixed freeze-out temperature 
T f (Wanajo et al. 2002, Wanajo 2007). Winds are charac- 
terized by a monotonically increasing velocity and a contin- 
uous temperature decrease when the radius goes to infin- 
ity, whereas breezes are obtained when a prescribed, non- 
vanishing value of the pressure and temperature is required 
to be reached at the outer boundary at some large radius 
(this is supposed to mimic the fact that the fast neutrino- 



driven wind is decelerated again as it collides and merges 
with the preceding, more slowly expanding material ejected 
in the early phase of the supernova blast). The approach 
taken by Wanajo et al. (2002) allowed them to causally 
disconnect the conditions at large radii from those at the 
neutron star and thus to choose the asymptotic tempera- 
ture independent of the mass-loss rate (in contrast to the 
situation for breezes and winds). The results from these dif- 
ferent types of outflow scenarios concerning combinations of 
entropy, Y Cl and expansion timescale necessary for strong 
r-processing in the ejecta turned out to be in qualitative 
agreement (compare, e.g., Fig. 11 in Thompson et al. 2001 
with Fig. 8 in Otsuki et al. 2000 and with Table 2 and Fig. 6 
in Sumiyoshi et al. 2000). 

Later Terasawa et al. (2002) announced to have found 
successful conditions even for proto-neutron stars with a 
more typical mass around 1.4 Mq and a radius of 10 km. 
Like previous studies (in particular Sumiyoshi et al. 2000, 
but also Otsuki et al. 2000) they considered breeze outflows, 
but different from the earlier investigations they chose a 
smaller value for the outer boundary pressure, which im- 
plied a lower asymptotic temperature. They argued that 
this is favorable for a strong r-proccssing up to the third 
peak because the lower final temperature is associated 
with a faster expansion and more rapid cooling. The more 
quickly decreasing temperature leads to a slowing down of 
charged-particle reactions and reduces the efficiency of a- 
particle recombination. Therefore it leads to less produc- 
tion of seed nuclei and a higher neutron-to-seed ratio. Such 
a sensitivity to the expansion timescale was also seen by 
Arnould et al. (2007), who performed systematic variations 
of the parameters in analytic wind and breeze solutions. 
Arnould et al. (2007) verified that wind solutions provide 
the more favorable conditions for strong r-processing than 
the slower breezes, because for given values of the entropy 
and mass loss rate the expansion timescale is directly cor- 
related with the asymptotic value of the temperature and 
thus of the total specific energy of the outflowing gas. These 
values are lowest in case of freely expanding winds. Arnould 
et al. (2007), however, also saw that the mass loss rate has a 
much more sensitive influence: breeze solutions with higher 
mass loss rates (and otherwise the same characterising pa- 
rameter values) make a faster expansion and allow for a 
stronger r-process despite having higher asymptotic values 
of the temperature. 

Although Arnould et al. (2007) confirmed the formation 
of the third r-process peak for combinations of entropy, 
expansion timescale, Y c , neutron star mass, and asymp- 
totic temperature in the ballpark of those considered by 
Terasawa et al. (2002), the results of the latter paper 
are nevertheless in contradiction to the earlier studies by 
Otsuki et al. (2000), Sumiyoshi et al. (2000), and Thompson 
et al. (2001): Terasawa et al. (2002) obtained a signifi- 
cantly higher entropy and shorter expansion timescale than 
Sumiyoshi et al. (2000) even for the same choice of outflow 
determining conditions, i.e., for the same individual neu- 
trino luminosity (L Vi = 10 51 ergs -1 ), the same mean neu- 
trino energies, and in particular the same outer boundary 
pressure (compare Table 1 in the Terasawa et al. work and 
Section 3.2 in Sumiyoshi et al.). Moreover, for all tested val- 
ues of the boundary pressure, Terasawa et al. (2002) found 
outflow properties for their 1.4M neutron star that were 
largely different from those plotted for breezes in Fig. 8 
of Otsuki et al. (2000) and for freely expanding winds in 
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Fig. 11 of Thompson et al. (2001). These differences seem 
to have been causal for the successful solar system like r- 
process reported by Terasawa et al. (2002), but the actual 
reason why the more favorable outflow behavior was ob- 
tained, remains unexplained- 

Recently, Arcones et al. (2007) performed new hydrody- 
namic simulations of neutrino-driven winds, systematically 
exploring the effects of the wind termination shock that 
forms when the supersonic wind collides with the slower 
earlier supernova ejecta (Janka & Miiller 1995a, b; Burrows 
et al. 2005; Buras et al. 2006, see also Tomas et al. 2004) 
and that decelerates the outflow abruptly. They found that 
the position of the reverse shock is strongly dependent on 
the evolution phase, progenitor structure, and explosion en- 
ergy of the supernova. Motivated by these studies, interest 
has recently turned to a closer exploration of the relevance 
of the late-time wind dynamics for r-process nucleosynthe- 
sis. On the one hand, Wanajo (2007) showed numerically 
that a solar-like r-process can also be produced in superson- 
ically expanding outflows whose temperature drops quickly 
to a few 10 8 K instead of asymptoting to a value around 
10 9 K as previously mostly assumed (Otsuki et al. 2000, 
Sumiyoshi et al. 2000, Wanajo et al. 2001, 2002). In such 
a low-temperature environment an (n, 7)-(7, n) equilibrium 
is never achieved during the nucleosynthesis of heavy r- 
process material, but neutron captures compete with (3- 
decays in the low-density matter, similar to what was dis- 
cussed by Blake & Schramm (1976). 

On the other hand, Kuroda et al. (2008) started to ex- 
plore for the first time systematically the consequences of 
the wind termination shock for the r-processing in the wind. 
Decelerating the outflow abruptly, the reverse shock does 
not only raise the entropy of the matter, but in particular it 
slows down the temperature and density decline that takes 
place in the subsequent expansion. Kuroda et al. (2008) 
pointed out that the change of the temperature behavior 
plays a decisive role in determining the r-process abun- 
dances, because the nucleosynthesis path depends strongly 
on the temperature during the r-process freeze-out phase. 
In contrast, the entropy jump does not seem to be im- 
portant because high entropies in the shocked outflows are 
reached only when the temperature is already well below 
2.5 x 10 9 K, i.e. not in the regime where the neutron-to-seed 
ratio is established before the onset of r-processing. 

In the present paper we also investigate the influence 
of the late-time outflow dynamics on the r-process nucle- 
osynthesis. To this end we consider outflow trajectories 
that consist of an initial homologous phase and a second, 
slower expansion stage. This describes the wind dynamics 
only schematically but at least some basic features that are 
found in detailed solutions of unshocked as well as shocked 
outflows can be reproduced. In contrast to the outflow so- 
lutions studied by Arnould et al. (2007), for example, our 
simple parametric ansatz for the temperature and density 
decrease in the ejected matter allows us to modify the early 
expansion behavior and that at late times independently, 
i.e. in an uncoupled way. This is similar to the approach 



1 Note that Thompson et al. (2001) used a sophisticated de- 
scription of the equation of state and thus explicitly accounted 
for the nonrelativistic character of electrons and positrons at low 
temperatures, an effect that Sumiyoshi et al. (2000) pointed out 
to be important for a reliable determination of the expansion 
timescale. 



taken by Wanajo et al. (2002) and closer to how, for ex- 
ample, the presence of a wind termination shock acts on 
the matter. Motivated by the insensitivity to the entropy 
jump seen by Kuroda et al. (2008), we ignore the disconti- 
nuity of the fluid variables at the location where their time 
dependence is assumed to change. 

In order to obtain r-processing with reasonable instead 
of extreme values of the wind entropy s and electron frac- 
tion Y" c , we consider sufficiently short exponential timescales 
for the first expansion stage, i.e., we consider timescales 
in the ballpark of those given for example by Otsuki et 
al. (2000) for outflows from a neutron star with a gravi- 
tational mass of 2Mq, a radius of 10km, and a neutrino 
luminosity of 10 52 ergs _1 (see their Table 1). We will vary 
the prescribed expansion behavior during the second phase 
in order to discuss the effects of the late-time temperature 
evolution on the r-processing from a more general point of 
view than done in previous works. The role of a strong or 
weak deceleration of the very fast, transsonic winds will be 
investigated by a constant or power-law behavior, respec- 
tively, of the temperature and density during the second 
stage. The former may be interpreted as the limiting case of 
a strong deceleration by a wind termination shock. The lat- 
ter may be considered as an approximation to the reduced 
acceleration that a freely expanding wind experiences at 
larger radii after its roughly homologous initial phase. 

Our paper is structured as follows. In Sect. 2 we will 
describe the dynamical model and the network used in our 
study, and will present some numerical tests we performed. 
In Sect. 3 we will describe our results and in Sect. 4 we 
finish with conclusions. 



2. Network, input data, and numerical modeling of 
the (a+r)-process 

We consider here the conditions for r-process nucleosynthe- 
sis in neutrino-driven outflows from the surface layers of hot 
nascent neutron stars, which have been the subject of many 
previous studies, where different aspects of the problem 
were discussed (see, e.g., Qian & Woosley 1996; Hoffman et 
al. 1997; Woosley et al. 1994; Cardall & Fuller 1997; Qian 
& Wasserburg 2000; Takahashi et al. 1994; Otsuki et al. 
2000; Sumiyoshi et al. 2000; Thompson et al. 2001; Kuroda 
et al. 2008, Arnould et al. 2007, and references therein). 
Our main goal here is to evaluate numerically the possi- 
bility of producing r-process elements in this environment 
in dependence of the late-time behavior of the outflowing 
gas. Within the framework of a very simple, purely ana- 
lytic model of the wind we want to determine the favorable 
combinations of entropy, Y e and dynamical timescale, and 
how sensitively these parameters influence the outcome. 
Though, of course, our results will not allow us to make 
a judgement about whether neutrino-driven winds are the 
long-sought site of r-process material or not, our paramet- 
ric approach has the advantage of reducing the dynamical 
aspects of the problem to an absolute minimum of ingredi- 
ents, giving one much freedom in the choice of the involved 
parameter values. 

2.1. Outflow behavior 

We represent the outflow behavior during the early and late 
expansion phases by different analytic functions. These de- 
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scribe qualitatively (but certainly not to very high accu- 
racy) on the one hand the wind acceleration through the 
sonic point, and on the other hand the evolution of the 
outflow during a second phase of either deceleration by a 
reverse shock or reduced acceleration. 

In the first stage of the expansion of a spherical mass 
shell we assume a homologous velocity-radius dependence, 
v oc r, corresponding to an exponential growth of the ra- 
dius, r(t) — i?i„i exp(i/rd yn ). Therefore, steady-state con- 
ditions (which imply r 2 pv = const) yield an exponential 
decline of the density and for an adiabatically expanding, 
radiation-dominated wind (i.e., p oc T 3 ) also an exponential 
decrease of the temperature: 



P(t) 
T 9 (t) 



Pini • cxp(-3t/r dyn ) 
exp(-t/r dyn ) . 



rplYll 

— 1 9 



(i) 

(2) 



Here Tg is the temperature normalized to 10 9 K and pi n i 
and Tg nl are the initial values of density and temperature 
at some small radius i?i n i. The dynamical timescale or ex- 
pansion timescale Td yn will be treated as a free parameter 
and can vary between 1ms and more than 100 ms. 

For reasons of simplicity, we will always assume that 
p is proportional to T 3 and that the gas entropy per nu- 
cleon (in units of Boltzmann's constant fee) is given by 
the relation s = 3.3AT^/p 5 with p 5 = p/(10 5 gem -3 ). One 
should, however, keep in mind that this is only a simpli- 
fying approximation, which is accurate only when the en- 
tropy of the gas is large enough (s <; 100 fcs per nucleon) 
and the temperature sufficiently high (Tg J> 5). In this case 
radiation-dominated conditions prevail, electrons are rela- 
tivistic, electron-positron pairs are abundant, and baryons 
contribute to the total entropy only at the level of a few 
percent. In reality, however, the pairs begin to disappear 
below Tg « 5 and the ratio T 3 /p increases. Assuming it 
to be constant also at low temperatures therefore leads to 
an overestimation of the density in the outflowing matter 
compared to a truely adiabatic evolution (see Witti et al. 
(1994) for a detailed discussion). 

The deceleration by the reverse shock or changing ac- 
celeration behavior are supposed to happen at a radius ro 
and time to when the velocity, density, and temperature 
reach the values vo, po, and To, respectively. As discussed 
in the introduction, we ignore here the discontinuous be- 
havior of the dynamical and thermodynamical variables at 
the shock and connect the early-time behavior continuously 
with the late-time behavior. This means that also the ra- 
tio T 3 / p (and therefore the quantity s that we consider 
as gas entropy) remains unchanged at the transition point 
and is taken to be constant during the following second 
stage of the expansion. For the latter we consider two cases 
with different limiting behavior for t^> to- In the first case 
we assume that the density and temperature asymptote to 
constant values, 



p(t) = Po, 

T(t) = T , for t > t . 



(3) 
(4) 



For steady-state conditions this implies that the radius and 
the velocity of a Lagrangian mass shell evolve at t > to 
according to 



r(t) 



ro 



1 + 3 ^ (t - 1 ) 
ro 



-1 1/3 



(5) 



v(t) 



Vo 



1 + 3 — (t- to) 
ro 



-2/3 
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and therefore r(t) oc i 1 / 3 and v(t) oc i~ 2 / 3 oc r~ 2 — » 
for t ^> t . We consider this choice of time-dependence as 
an approximative representation of the strong deceleration 
experienced by the wind in passing through a termination 
shock. Assuming constant density and temperature on the 
typical timescale of r-processing (several hundred millisec- 
onds) is in reasonably good agreement with the slow evolu- 
tion of the shocked outflow found in detailed hydrodynamic 
models (see Figs. 6 and 8 in Arcones et al. 2007). 

In the second investigated case the density and tempera- 
ture are still assumed to decline at late times, but much less 
steeply than during the exponential first expansion phase: 



Pit) = poi- 



nt) = To 



-2/3 



, for t > to ■ 



(7) 

(8) 



For steady-state conditions this corresponds to 
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v(t) = v 
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To 



(9) 
(10) 



which yields r(t) oc t and v(t) — v^ 3 (ro/to) 2 ^ 3 — const > 
for t 3> to- A posteriori this asymptotic behavior justifies 
our choice of the power-law time-dependence for p and T in 
Eqs. ((TJ) and ©. Note that v(t — > oo) < i>o, i.e. deceleration 
(and not just a slow-down of the acceleration) happens only 
at late times if to > Td yn - Values of the parameters used for 
some of the considered model cases are listed in Table [TJ 

The slower decline of the density and temperature dur- 
ing the power-law phase may be considered as a simplified 
representation of a reduced wind acceleration occurring af- 
ter an approximately homologous (i.e., v oc r) early expan- 
sion. Neutrino-driven winds exhibit such a nearly homolo- 
gous evolution and thus nearly exponential increase of the 
radius and velocity with time only up to some distance, but 
then the velocity continues to grow less rapidly even in the 
absence (or before) a possible strong deceleration by a ter- 
mination shock. The basic properties of this behavior are 
captured by our approach. 

2.2. Calculating the nucleosynthesis 

At the onset of its expansion, the ejected matter in the 
neutrino-driven wind is very hot and composed of free neu- 
trons and protons. With ongoing cooling the nuclear statis- 
tical equilibrium (NSE) shifts towards an increasing mass 
fraction of alpha particles until finally the recombination 
to heavy nuclei sets in. Provided the conditions of entropy, 
electron fraction, and expansion timescale are suitable, the 
production of elements between the first and second peaks 
of the abundance curve through charged particle reactions 
and neutron captures may occur. For sufficiently extreme 
conditions even third-peak elements and (A ~ 196) may 
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be assembled. Typically starting our nucleosynthesis cal- 
culations at temperatures Tg ~ 6 and densities p ~ 10 5 - 
10 6 g cm -3 (the initial density is chosen according to a spec- 
ified value for the conserved gas entropy), the element for- 
mation in our model runs first proceeds mainly through 
charged-particle reactions. As the temperature and density 
decrease, the importance of (a,X)- and inverse reactions di- 
minishes and the nuclear flow begins to be driven by a dy- 
namical r-process, provided a sufficient number of free neu- 
trons is still present. The production of a small mass frac- 
tion of heavy nuclei with Z > 26 before the onset of rapid 
neutron captures significantly reduces the requirements for 
the neutron source (a smaller number of free neutrons is 
required) . 

The time-dependent concentrations of nuclear species, 
Y(A,Z), during the r-processing as determined by reac- 
tions with neutrons, beta-decays, and fission processes are 
described by the nucleosynthesis network developed by 
Blinnikov & Panov (1996) and Nadyozhin et al. (1998). 
This network was here extended by charged-particle reac- 
tions and a larger number of nuclear reactions and fission 
processes so that it is possible to handle both the initial a- 
process and the subsequent r-process with the same code. 

The number of nuclei and reaction equations included 
in calculations depends on boundary conditions and the 
employed nuclear mass model and can be as large as about 
4300. We considered the region of nuclei with Z ranging 
from 3 to 100. The minimum and maximum atomic mass 
values for each chemical element were determined by the 
proton and neutron drip lines. 

The reaction rates entering the system of differential 
equations differ by tens of orders of magnitude. Thus, the 
system of equations for nuclear kinetics to be solved is a 
classical example of a stiff system of ordinary differential 
equations. We used one of the most effective methods to 
integrate such a stiff system of equations, Gear's method 
(Gear, 1971). The description of the complete package of 
solver routines and its applications to the r-process calcu- 
lations can be found in Nadyozhin et al. (1998). 

Nuclear mass values as predicted by FRDM (Moller et 
al. 1995) were used, the beta-decay rates were calculated in 
the framework of the QRPA-model (Kratz et al. 1993), and 
the reaction rates with neutrons were described according 
to the calculations of Cowan et al. (1991) and of Rauscher 
& Thielemann (2000). 

In the nucleosynthesis studies presented here, the triple 
a and aan reactions of helium burning, 3a — > 12 C and 
a + a + n — > 9 Be, respectively, along with their inverse re- 
actions, were included. The rates for both processes were 
taken from the REACLIB library of Thielemann et al. 



(1987) and those for the reactions of heavier nuclei with pro- 
tons and a-particles from Rauscher & Thielemann (2000). 

2.3. Code tests and comparisons 

Code test s for a number of different cases were per- 
formed by Panov et al. (2001a)| and for a number of ex- 
plored cases r-process calculations with the same rates 
ga ve practically the same iso topic yields as in the paper 



of Freiburghaus et al. (1999) Network calculations of the 



a-process with the present code were compared with the 
results shown by Witti et al. (11994p and yield rather com- 
patible abundance distributions at the beginning of the r- 
process, in spite of differences in the employed nuclear re- 
action rates and the nuclear mass model. 

In this context we would like to note that the transition 
from the a-process to the r-process can only be done cor- 
rectly on the basis of the same mathematical model. In our 
calculations we use the same code for both parts of the nu- 
cleosynthetic reaction sequence, without artificial devision 
into a- and r-process steps. 

We compared the seed production as calculated with our 
full network to results from nuclear statistical equilibrium 
(NSE) as obtained by Nadyozhin & Yudin (2004), Blinnikov 
et al. (to be published), and with the NSE-code used by 
Witti et al. ( 1994). We observed rather good agreement at 
temperatures Tg « 5-6 and a density around 5 x 10 6 g cm~ 3 . 
At such temperatures the a-peak is already formed and the 
mass fraction of seed nuclei (i.e., heavy elements) is still 
small, the most abundant nuclei in both cases being 50 Ti 
and 54 Cr. Small discrepancies between the different calcula- 
tions emerge mostly from the use of different mass formulas 
and can lead to smaller differences in the subsequent a- and 
r-process. They may be important for exact predictions of 
the abundance flow during the r-processing, but they should 
not affect the basic results of our investigation. 

We also compared our results for the a-process with the 
calculations made by Witti et al. ([1994), using conditions 
similar to those of the HT-model of Witti et al. £H)93 h i.e., 
s = 390, Y c = 0.455. For an exponential expansion with 
a timescale of Td yn = 62 ms, which is the same as in the 
HT-model, the number density function of tem- 

perature (Fig. [TJ left) is clearly different from what was 
obtained by Takahashi et al. I|1994[) (see Fig. 2 therein and 
the dashed line in the left panel of Fig. [I}. When we in- 
crease the expansion timescale to Td yn = 162 ms, our result 
becomes close to that of Takahashi et al. (1994]). The reason 
for this discrepancy is the fact that our calculations assume 
homologous expansion with an exponential decrease of den- 
sity and temperature, whereas the temperature and density 
in the hydrodynamic model considered by Takahashi et al. 
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Fig. 1. Tg-dependence of the free neutron density n n (in particles per cm -3 ) for s — 390 (measured in k-Q per nucleon), 
Y c = 0.46 (left) and s = 145, Y c = 0.42 (right). The different curves correspond to different exponential expansion 
timescales Td yn as given in the list. The dashed line in the left plot shows the result for conditions close to those in Fig. 7 
of Witti et al. (|1994p . In the case of the HT- model considered by Witti et al. (1994), the temperature asymptotes to a 
value of Tg = 1.5. The density also decreases so slowly that a high neutron number density is supported long enough 
for the creation of a strong third r-process peak. In contrast, during the free expansion of very high entropy material 
with an exponential timescale of 0.62 ms, the matter dilutes too fast to produce large amounts of heavy elements, and 
no prominent third r-process abundance peak can be formed. In spite of a similar decline of the neutron density with 
temperature for an expansion timescale of 162 ms, the third abundance peak can also not be assembled in this case, 
because the slow expansion leads to abundant seed production in a strong a-process and too small a neutron-to-seed 
ratio. 



asymptote to nonvanishing values, i.e., the temperature ap- 
proaches Tg — 1.5 during the later expansion. Although 
r-process elements are formed, neither of our exponential 
runs leads to the build-up of a strong third r-process peak, 
despite the high entropy. While in the case with the short 
expansion timescale (rd yn = 62 ms) the rapid dilution of 
the matter prohibits the efficient formation of seed nuclei, 
the larger expansion timescale (Td yn = 162 ms) leads to a 
strong a-process with a lot of seed production and there- 
fore a neutron-to-seed ratio that is too small for the gener- 
ation of very heavy r-process nuclei. In contrast, Takahashi 
et al. (|1994p obtained a prominant third abundance peak. 
In their calculation a sufficiently fast initial expansion pre- 
vents the formation of too much seed material before the 
r-process starts, and therefore the neutron-to-seed ratio re- 
mains high. The asymptoting temperature and density on 
the other hand support a high neutron number density for 
such a long time that neutron capture reactions can assem- 
ble nuclei also in the third peak. 

For lower entropy, s — 145, and higher neutron ex- 
cess, Ye = 0.42, (Fig. |TJ right) the number of free neutrons 
can also remain large enough to allow for the onset of re- 
processing when the a-process freezes out at a temperature 
of Tg ~ 2. In this case, however, much shorter dynamical 
timescales are needed. Figure [5] displays the decreasing seed 
formation and increasing neutron-to-seed ratio for smaller 
expansion timescales ra yn . 

How many neutrons per seed nucleus are needed to form 
the third abundance peak (A ~ 196)? Usually the number 
mentioned in this context is not less than 150. But this esti- 
mate is based on the simple calculation how many neutrons 
a single 56 Fe nucleus must capture to finally, after a chain 
of beta-decays, end as nucleus in the platinum peak. In the 



actual r-process, only a fraction of the nuclei that are ini- 
tially formed by the alpha-process — we call these nuclei 
the "seed" — will be ultimately transformed to elements 
and isotopes with the highest mass numbers. 

Taking into account that the observed ratio of the third 
to second abundance peak is Yig6/Fi3o ~ 0.2, we estimate 
that the neutron number needed to obtain a strong third 
peak is of the order of (196 - 60) x 0.2 w 30. The actual 
value might even be a bit lower, because the mean atomic 
number at the freeze-out time of the a-process can be about 
80 or more instead of 60 (an exact estimate, however, also 
depends on the conditions in which the r-processing takes 
place and the corresponding speed of the abundance flow). 
Therefore a solar-like formation of the platinum peak can 
be expected if one has a neutron-to-seed ratio of around 30 
after the freeze- out of charged-particle reactions. 

We point out, however, that this value does not ensure 
that the solar abundance distribution in the whole mass- 
number region between the second and third abundance 
peaks is reproduced well. This goal is not attempted in our 
study here, and it is well known that a superposition of at 
least two components with different physical conditions in- 
stead of a single ejecta trajectory (as considered by us) is 
needed for matching the solar r-process pattern (see, e.g., 
Kratz et al. 1993, Goriely & Arnould 1996, Wanajo et al. 
2004, Kuroda et al. 2008). In our fast-expansion cases we 
regard the described criterion as a reasonably good indi- 
cator for outflow conditions that allow for the appearance 
of a strong third peak during the r-processing that follows 
after the freeze-out of charged-particle reactions. 

To judge about the possibility of an r-process for dif- 
ferent choices of the parameters of our dynamical model 
(s, Y c , Td yn ), we first consider the seed formation and the 
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Fig. 2. Tg-dependence of the neutron-to-seed ratio Y n /Y s (left) and of the "seed" number fraction Y s (right) for an 
entropy of s — 145 and Y c = 0.42. The different curves correspond to different exponential expansion timescales, 
Tdyn = 5, 10, 20, 30, 50 ms (the same cases as in the right panel of Fig. [1]). At the beginning of r-processing (Tg » 2), the 
free neutron number is sufficiently large to form heavy elements up to the second or even third peak if the dynamical 
time is less than 10 ms for the considered entropy and electron fraction. For larger values of Td yn , the free neutrons are 
exhausted right at the beginning of the r-processing and the number of seed nuclei is too high. 



corresponding time evolution of the neutron-to-seed ratio 
in the homologous expansion phase. In Fig. [T] we have 
seen agreement of the free neutron density as a function 
of time between our calculations and the prior ones by 
Takahashi et al. (1994) for a suitable choice of the expansion 
timescale Td yn (although the heavy-element nucleosynthe- 
sis was considerably different as discussed above). Figure^ 
shows the neutron-to-seed ratios and the seed abundances, 
Y s — J2z>2^z, versus temperature for the same cases as 
displayed in the right panel of Fig. [U Note that we consider 
all nuclei with Z > 2 as seeds, which is a different definition 
than used by Terasawa et al. (|2001j) . who restriced seed nu- 
clei to the more narrow range of 70 < A < 120 and Z > 26. 
The latter mass range is appropriate when the formation 
of the second abundance peak in an incomplete r-process 
is discussed. Here, however, we explore the possibility of 
third-peak production and consider a combined a- and r- 
process. In this case the atomic mass number range should 
not be constrained. Because of the very short expansion 
timescale and the treatment of both the a-process and the 
r-process within the framework of the same mathematical 
model and computer code, a distinction of seed nuclei and 
r-processed material becomes artificial and one has to care- 
fully judge when this information is measured and what 
it means for the evolution of the mass number with time. 
Often it refers to older calculations in which the a- and r- 
process were computed in two independent steps and with 
different codes. In this context we note that the amount 
of heavy nuclei in our calculations is approximately twice 
as big as found by Terasawa et al. (2002). This difference 
is mostly explained by the different model parameters (in 
particular, Terasawa et al. investigated outflows with higher 
entropies), but to some extent may also be a consequence 
of different nuclear rates and included reactions. 

With the definition of seed nuclei adopted by us, the 
seed abundance becomes constant after the freeze-out of 
the alpha-process and will also be unaffected by the transi- 
tion to a second stage with modified expansion behavior for 



the cases considered by us. In contrast, replacing the free 
(homologous) expansion by a slower second stage of outflow 
behavior will have an impact on the density of free neutrons 
as a function of time. The free neutron density determines 
the r-process path and thus the formation of the heaviest 
elements in the third abundance peak. 

From Fig. fj] (left panel) we see that the expansion 
timescale for third peak formation to happen should be 
less than about 10 ms. In this case the neutron-to seed ratio 
Y n /Y s as the decisive macroscopic factor for the platinum 
peak formation reaches the interesting values mentioned 
above. Naturally, this ratio depends not only on the number 
of free neutrons, but also on the seed abundance. Figure [5] 
(right panel) shows that the seed production varies strongly 
with the expansion timescale: if the dynamical time is large, 
a lot of seed is assembled and all neutrons will be cap- 
tured faster than the platinum peak begins to appear. With 
shorter expansion timescale the seed production drops and 
at the same time the density of free neutrons increases (see 
right panel of Fig. [TJ . The consequence of both trends with 
reduced Td yn is a strong growth of the Y n /Y s ratio that is 
present at the beginning of the r-process. 

However, as already discussed above in the context of 
Fig. [2 too rapid expansion can also be disadvantageous for 
a strong r-process. When the dynamical timescale becomes 
very short, the expansion and gas dilution proceed faster 
than the rate of recombination of alpha particles to seed 
nuclei and subsequently the rate of neutron captures. For 
the conditions considered in Fig. [2] this happens when the 
homologous expansion and density decrease are too rapid 
to allow for the neutrons to be assembled into a strong 3rd 
peak, despite a high neutron-to-seed ratio and a sufficiently 
small seed concentration Y s below a critical upper value 
of Y^ 1 ' ~ 3 x 10 -3 after the freeze-out of charged-particlc 
processes. 

Of course, as discussed in detail in many previous 
works, besides the expansion timescale, the entropy and 
the neutron excess (or Y c ), have a sensitive influence on 
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Fig. 3. Time evolution of the temperature for the dif- 
ferent considered cases with exponential first expansion 
phase (starting at Tg ~ 6 and declining with a dynami- 
cal timescale of Td yn = 5 ms) and constant later phase. 



the strength of the r- processing, i.e., on the question how 
many nuclei with mass numbers above the second abun- 
dance peak and in particular near the third peak can be 
formed. In any case, for rapid expansion (small Td yn ) the 
seed production is reduced and the neutron-to-seed ratio 
becomes more favorable for a strong r-process. The heavier 
the seed nuclei at the end of the a-process are, the lower can 
the required number of free neutrons be. This was already 
discussed by, e.g., Panov & Chechetkin (2002), who showed 
that fairly low neutron/seed ratios are already sufficient 
when the seed material at the freeze-out time of charged- 
particle reactions consists mostly of nuclei in the second 
abundance peak. 



3. Asymptotic behavior of temperature and density 
and formation of the platinum peak 

Our combined (a+r)-code was applied to nucleosynthe- 
sis calculations assuming the two-stage expansion behav- 
ior described above, with a second phase of either constant 
or slowly decreasing temperature and density following a 
first phase of rapid, exponential expansion. Our calcula- 
tions were started at a temperature of Tg = 6, assuming 
NSE at this point. For conditions similar to those given 
in Fig. 2 of Terasawa et al. (2001) with a short expansion 
timescale of Td ya = 5 ms, we obtained basically the same 
results for the development of the neutron number density 
and average neutron separation energy until about half a 
second, and found gross agreement of the structure of the 
abundance distribution. Smaller discrepancies might be at- 
tributed to differences in the initial composition and nuclear 
reaction rates. The number fraction of seed nuclei reported 
by Terasawa et al. (2001), Y s ~ 0.001, is a bit low com- 
pared to our results, but this is probably mainly caused by 
a different definition of "seeds" . 

In contrast, we were not able to confirm the formation 
of the third abundance peak for conditions similar to those 
considered by Terasawa et al. (2002) with an assumed ex- 
pansion timescale of Td yn = 25 ms. Besides not providing 
exact information about the initial density and the kind of 



network used for the alpha- and r-process calculations (a 
full network or off-line calculations for the r-processing?), 
some of the results and explanations are hard to reproduce 
in detail. For example, the calculations done by us with dif- 
ferent NSE codes (including that of Nadyozhin and Yudin 
2004) and also the simulations by Witti et al. (1994) show 
that a-particles reassemble mainly down to temperatures 
around Tg = 5-6, somewhat dependent on the density, but 
not at temperatures as low as Tg = 4 (see Fig. 1 in Terasawa 
et al. 2002). Moreover, they argued that a lower value of 
the asymptotic temperature is favorable for a successful 
r-processing by allowing for a higher neutron-to-seed ra- 
tio, because such a lower "outer boundary" temperature 
reduces charged-particle reactions and thus the produc- 
tion of seed nuclei. This suggested influence of the chosen 
asymptotic temperatures in the range between Tg = 0.4 
and Tg = 1.3, however, is implausible, because charged- 
particle reactions become inefficient already at tempera- 
tures of Tg w 2, when the thermal energies of protons and 
a-particles become too low for enabling these particles to 
overcome the nuclear Coulomb barriers. 

We will therefore not further attempt to compare our 
calculations with those done by Terasawa et al. (2002). 
Instead, we will in the following present our results for the 
formation of the third abundance peak in dependence of 
entropy, electron fraction, and dynamical timescale during 
the first, exponential expansion phase. In particular we will 
also study variations with the asymptotic values of tem- 
perature and density during the second, slower phase of 
expansion. This will help us in the analysis of the physical 
effects that can explain our results. Our goal is to develop 
a deeper understanding of the influence of the late-time be- 
havior of neutrino-driven winds on the possibility of strong 
r-processing in such an environment. 

To this end we carried out a set of calculations for ex- 
ponential timescales in the range from 1.0 to 25 ms and for 
four values of the wind entropy (in units of Boltzmann's 
constant per nucleon), s = 105, 145, 170, and 200. The ini- 
tial neutron to proton ratio was determined by an electron 
fraction of usually Y c — 0.42; some runs were performed 
with a value of Y e — 0.46. 

Our special attention was on the second stage of slower 
expansion, which was chosen to either proceed with con- 
stant temperature and density (Eqs. O IHl) or gradually de- 
creasing temperature and density (Eqs. [Tl- llOp . In the for- 
mer case we defined the asymptotic temperature as Tg(t > 
tg) = Tg — const (Sect. l3"7Tj) . in the latter case the temper- 
ature was assumed to follow a power-law time-dependence 
according to T 9 f (t) = T 9 f (i ) x (tg/t) 2 / 3 for t > tg (Sect. [3121). 
The range of temperatures Tg(tg), where the second expan- 
sion phase began, had a broad overlap with the "boundary 
temperatures" considered by Terasawa et al. (2001, 2002) 
and Wanajo et al. (2002). 

Our standard calculations were performed with an en- 
tropy of s = 145. Values of less than s ~ 100 turned out 
to lead only to the production of the second abundance 
peak around A ~ 130, but r-processing up to the platinum 
peak was not possible when the other characteristic wind 
parameters were varied within the limits mentioned above. 

3.1. Constant asymptotic temperature and density 

In this section we consider an exponential first expansion 
phase that is superseded at evolution time tg by a second, 



Panov & Janka: Proto-neutron star winds and r-process nucleosynthesis 



9 



le+28 



le+26 



le+24 



le+22 - — 



le+20 





Fig. 4. Time dependence of the neutron number density, n n (t) (left, in particles per cm~ 3 ), and the final abundance 
distributions, Ya, resulting from our network calculations for s — 145, Y c = 0.42, a short dynamical timescale of 
T dyn = 5ms, and different values of the constant asymptotic temperature Tg(t > to) = Tg(to) during the second stage of 
the ejecta expansion. The corresponding values of Tg are listed in the inset and the temperature evolution of function of 
time is shown in Fig. [3] 



slow expansion phase whose asymptotic velocity at large 
radii goes to zero. In this case the density and temperature 
during the second phase adopt constant values (Eqs. [3H5J). 
We explore various choices of the asymptotic temperature 
Tl(t > t Q ) = T 9 ° = const between 0.1 and 1.4 (see Fig.©. 
For an entropy of s = 145, for which results are displayed in 
Figs. 0]-[6l these temperatures correspond to constant den- 
sities po between about lgcm -3 and 10 4 gcm -3 . For these 
values the neutron number density at time to shows differ- 
ences by more than three orders of magnitude (see Fig.[4|. 
Although the lowest assumed value of Tg = 0.1 might ap- 
pear extreme, because it requires very rapid exponential 
expansion for a longer time with fairly high velocities at the 
end of this phase (see Tabic [TJ) , the wide range of asymp- 
totic temperatures allows us to better understand the de- 
pendence of the nucleosynthesis on the late-time expansion 
behavior of the outflowing matter. 

When the asymptotic temperature is as high as Tg > 
1.4, the third peak hardly develops (see the abundance dis- 
tribution at the end of our calculations in the right panel 
of Fig. [3]). Only when the asymptotic temperature is re- 
duced from this value to smaller numbers, the platinum 
peak grows in strength, because the free neutrons are ex- 
hausted more rapidly and the r-processing proceeds faster 
towards the high mass-number region. The latter fact can 
be seen by comparing the decrease of the neutron density 
with time for Tg = 1.4, 1, and 0.4 in the left panel of 
Fig. IH and it is also visible from the neutron-to-seed ra- 
tios Y n /Y s as functions of time in the left panel of Fig. [5] 
When the asymptotic temperature is lowered to less than 
Tg w 0.3, the free neutron density n n during the slow- 
expansion phase becomes so low that the formation of the 
third abundance peak under such conditions slows down 
considerably. Correspondingly, n n as well as Y n /Y s decrease 
less quickly (see again the left panels of Figs. Q] and [5]) and 
the height of the third peak at the time of neutron exhaus- 
tion becomes clearly lower (Fig. 21 right panel). 

This inversion of the third peak formation with decreas- 
ing value of the asymptotic temperature Tg can be clearly 



seen from the time evolution of the average of the yields in 
the platinum peak relative to those in the N = 82 peak, 
(Y(A 3 )/Y(A 2 )), in the right panel of Fig. H The displayed 
quantity is defined by summing up the produced yields for 
five mass numbers around A — 130 and A = 196 and then 
computing the ratiq^l 

The left panel of Fig. [6] shows the sensitivity of the fi- 
nal value of this ratio to the chosen temperature Tg(t > 
to) = Tg. While the growing strength of the platinum 
peak with lower asymptotic temperature agrees with the 
finding by Terasawa et al. (2002), their conclusion of a 
more efficient r-processing for lower "outer boundary pres- 
sure/temperature" obviously does not hold any more when 
this temperature is less than T 9 0.3, which is below the 
range explored by them. 

We note here that often the change of the average 
atomic number is used for following the assembling of heavy 
elements beyond the A » 130 abundance peak. We found 
this quantity to be less suitable for this purpose than our 
height ratio (Y (A3) / Y (A 2 )) , in particular when the second 
peak is much stronger than the third. A growth of the lat- 
ter by one order of magnitude can mean a change of (A) 
by just one or a few units, whereas (Y (A3) / Y (A 2 )) varies 
sensitively and thus serves well as a tracer of changes of 
the abundance distribution in the high-mass number region. 
The numbers listed in Table[TJ for example, confirm the use- 
fulness, in fact the superiority, of the ratio (Y (A3) / Y (A 2 )) 
instead of the mean mass number, (A), as an indicator of 
the third peak formation at least in the cases with very 
fast outflow expansion considered here (of course, again 
we do not claim that a single parameter value or a sin- 
gle outflow trajectory is sufficient for getting a good match 
of the whole r-process abundance distribution between the 

2 A similar peak ratio was considered in a recent paper by 
Beun et al. (2008). 
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Fig. 5. Time dependence of the neutron-to-seed ratios, Y n /Y s (left), and of the height of the third abundance peak relative 
to the second, (Y (A3) / Y (A2)) (right, in percent) for the same conditions as in Fig. |4] (s = 145, Y e = 0.42, Td yn = 5 ms, 
Tg = const as given in the inset). The inversion of the third peak formation with decreasing asymptotic temperature can 
be clearly seen: when Tg is lowered from 1.4 to 0.4, the ratio of the third peak to the second increases, and when Tg is 
reduced even further, {Y (A3) / Y (A2)) begins to drop again. 




Fig. 6. Left: The changing strength of the third abundance peak, measured in terms of the ratio (Y (A3) / Y (A2)) , as 
function of the asymptotic temperature Tg(t > tg) = Tg for the conditions considered in Fig. |4] (s — 145, Y c = 0.42, 
Tdyn = 5 ms, Tg = const). The relative height of the third abundance peak depends strongly on the value of Tg because of 
a sensitive influence of this temperature on the r-process path. This can be seen in the distribution of cadmium isotopes at 
an evolution time of t = 0.1 s for three different cases of the asymptotic temperature (right; the neutron number densities 
in the inset are those at t = 0.1s in Fig. [4|. When Tg is lowered from 1.4 to 0.4, the free neutron density and thus 
the neutron capture rate in the slow, second phase of the expansion decreases less strongly than the (7,n)-rates, which 
are extremely sensitive to the temperature. As a consequence, the r-process path is shifted towards the neutron-drip 
line (compare the isotope distribution given by the green line) where the /3-decay rates are higher and the r-processing 
proceeds so rapidly that there is time to form a strong third abundance peak. For even lower asymptotic temperature 
the (7,n)-rates do not play an important role and the strength of the r-processing is determined by the free neutron 
density and the /3-decay rates. Since the r-process path returns again to a location closer to the valley of stability (see 
the isotope distribution of the blue line), the /3-decay rates are lower and the nucleosynthesis slows down. Therefore a 
strong third peak cannot be formed. 



second and third peaks). While (Y (A3) / Y (A2)) exhibits a 
rapid growth during the build-up of the third abundance 
peak and changes by many orders of magnitude between 
cases with weak or strong third peak, the average mass 
number remains within the relatively narrow interval of 
118 < (A) < 145 (see Table [J). 



The variation of the strength of the third peak with 
different asymptotic temperatures can be understood from 
the sensitivity of the neutron capture rates and nuclear 
photodisintegration rates to the neutron number density 
and temperature, respectively, and by the competition of 
these rates. This competition determines the location of 
the r-process path and thus the speed of the nucleosynthe- 
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Fig. 7. Time evolution of the temperature for the dif- 
ferent considered cases with exponential first expansion 
phase (starting at Tg ~ 6 and declining with a dynami- 
cal timescale of Td yn = 5 ms) and power-law later phase. 



sis, which is defined by the /3-decay rates. When the tem- 
perature during the second, slow expansion phase is large 
(Tg J> 1.0), the (7,n)-reactions are very fast and the r- 
process path lies closer to the stability region than in case 
of smaller Tg . The r-process flow beyond the second peak 
is then rather weak because of the low beta-decay rates 
and a correspondingly slow progression of the nuclear flow. 
When the asymptotic temperature is reduced to Tg * 0.2- 
0.7, the (7,n)-rates decrease and the r-process path moves 
towards the neutron-drip line where the /3-decay rates are 
higher. Therefore the r-processing proceeds faster beyond 
the second peak, leading to a more rapid drop of the free 
neutron density and a more efficient third peak production. 
When the asymptotic temperature is lowered to less than 
T° « 0.2 during the late expansion phase, the neutron den- 
sities are very low so that the r-process path returns to 
a location closer to the /3-stable region. In this case the 
/3-decays again become slower and therefore the r-process 
nucleosynthesis decelerates and the third peak builds up to 
a smaller height. 

These movements of the r-process path are a conse- 
quence of the different influence of a change of the asymp- 
totic temperature on the (7,n)- and the neutron captures 
rates. For asymptotic temperatures in the interval 0.2 <, 
Tg 1.0 the (7,n)-rates are very sensitive to temperature 
variations. In contrast, the neutron capture rates change 
with the corresponding variations of the density and free 
neutron density less strongly. A reduced temperature there- 
fore decreases the photodisintegration rates significantly, 
whereas the neutron captures remain fast despite the lower 
neutron density. When Tg ^ 0.2, the (7,n)-rates become 
less relevant, but neutron captures still compete with (3- 
decays. For such low temperatures and thus low neutron 
densities, the neutron capture rates are too small to drive 
the r-process path far away from the valley of stability. 

The shift of the r-process path that is caused by differ- 
ent asymptotic temperatures (and connected parameters) 
manifests itself in different isotopic profiles of the elements 
formed by the r-processing (i.e. in different yield of the 



isotopes of an element). This can be seen for the case of 
cadmium in the right panel of Fig. [5] 

3.2. Power-law time-dependence of the asymptotic 
temperature and density 

In this section we consider the case that the initial expo- 
nential phase is superseded at time t = to by a slow ex- 
pansion phase in which the temperature and density de- 
cay according to power laws, i.e., T|(i) oc Tg(t )/t 2 ^ 3 and 
p f \t) oc p f (io)/i 2 for t > t (Eqs. I7HTD|1. The temperature 
evolution as function of time for the considered cases is 
displayed in Fig. [7J 

For our standard set of wind parameters, s = 145 (in 
units of Boltzmann's constant per nucleon), Y c — 0.42, and 
Tdyn = 5 ms already used in Sect. 13. ll Fig.[8]shows the time 
evolution of the free neutron density (left panel) and the fi- 
nal abundance distributions for four different values of the 
transition temperature Tg(to) (right panel). The formation 
of the third abundance peak turns out to be fairly insensi- 
tive to variations of Tg(to) between about 0.4 and 1.4. For 
all transition temperatures in this interval the temperature 
evolution is very similar (Fig. [7J and a prominent third 
peak appears. 

Since the temperature in the post-exponential phase 
drops rapidly, the strongly temperature-dependent (7,n)- 
reactions become unimportant very soon, while the con- 
sidered entropy allows a high free neutron density (n n > 
10 22 cm -3 ) to be present still for a long time. At such con- 
ditions the r-process path moves very close to the neutron- 
drip line and returns to the classical r-process path during 
free neutron exhaustion (Panov 2003). For the considered 
conditions its location shifts significantly only when the free 
neutron density changes by 2-3 orders of magnitude. This 
explains the relative robustness of the abundance yields to 
variations of Tg(to) around 1.0. 

The left panel of Fig. [5] displays the evolution of the 
neutron-to-seed ratios Y n /Y scc d that correspond to the neu- 
tron number densities of Fig.JS] With an entropy of s = 145 
and an exponential timescale Td yn = 5 ms the values of 
Y n /Y see d at the time when the r-process has made the 
second abundance peak (at t a to) are around 20 for 
Tg(t ) = 0.4-1.4. As discussed in Sect. [U this is a slightly 
insufficient number of neutrons per seed nucleus to create a 
third peak with exactly the observed yields. The small un- 
derabundance of the third peak can also be seen in the 
right panel of Fig. [5J where the height of the platinum 
peak relative to the tellurium peak is given as function of 
time. Different from the case of constant asymptotic tem- 
perature (see Fig. O, {Y(A 3 )/Y(A 2 )) is nearly the same 
for Tg(io) = 1-0 and 1.4, and is only slightly reduced for 
Tg(to) = 0.4. For even lower transition temperatures the 
height of the third peak drops steeply. This is evident from 
the blue line in the left panel of Fig. [TO] and is a consequence 
of a faster decrease of the free neutron density (compare 
the left panels in Figs. Hand [8]), which for T 9 l (t Q ) < 0.4 
falls quickly to a value of n n ~ 10 22 cm -3 , below which the 
equilibrium of n-captures and /3-decays is shifted to a re- 
gion of slow beta-decay rates, leading to reduction of the 
third peak. 

The left panel of Fig. [TO] provides a direct comparison 
of the variation of (Y (A3) / Y (A2)) with the transition tem- 
perature between the exponential first expansion phase to 
the second phase of either constant temperature or power- 
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Fig. 8. Same as Fig. [4] (with s = 145, Y c = 0.42, and an exponential timescale Td yn = 5 ms), but for a power-law decay 
of the temperature during the second, slow expansion phase: T|(t) oc Tg(to)/t 2 / 3 for t > to (see Fig. [7|). The left plot 
shows the time dependence of the neutron number density, n n (t) (in particles per cm -3 ), the right plot the corresponding 
final abundance distribution for different values of the transition temperature Tg(to) as given in the inset. Different from 
the case with constant asymptotic temperature in Fig. [4j the temperature and thus the photodisintcgration rates drop 
rapidly also in the slow expansion phase. Neutron captures are therefore able to produce a prominent third abundance 
peak even for the largest considered value of the transition temperature, T|(to) = 1-4. The height of the third peak 
relative to the second is very similar in the whole range of temperatures Tg(to) between 0.4 and 1.4, and the inversion 
behavior of (Y(A 3 )/Y(A 2 )) with T 9 f (i ) is absent. 




Fig. 9. Same as Fig. [5] (with s = 145, Y c = 0.42, and an exponential timescale Td yn = 5 ms), but for a power-law decay 
of the temperature during the second, slow expansion phase: Tg(i) oc Tg(tg) /t 2 ^ 3 for t > to (as considered in Fig. [5]). The 
left plot shows the time-dependent neutron-to-seed ratio Y n /Y a , the right plot the height of the third abundance peak 
relative to the second ((Y (A3) / Y (A 2 )) , in percent) as functions of time. The weak dependence of the r-processing on 
variations of the transition temperature Tg(tg) between 0.4 and 1.4 as visible in the right panel of Fig. [5] is also evident 
from these two plots. 



law temperature decrease. Above a transition temperature 
Tg(to) ~ 0.4 the behavior in both cases is dramatically 
different. While in the case of constant asymptotic temper- 
ature the relative height of the third peak decreases with 
higher values of Tg (see Fig. [6] and Sect. 13. 1|) . a power-law 
decay of the temperature in the second expansion phase 
leads to a prominent platinum peak for all values of Tg(to) 
between 0.5 and 1.4. The reason is again the large sensitiv- 
ity of the photodisintegration reactions to the late-time be- 
havior of the temperature. Because of the power-law decline 



the temperature drops within milliseconds to values where 
a (n,7)-(7,n) equilibrium is no longer possible but is re- 
placed by a quasi-equilibrium between (n,7)-reactions and 
/?-decays. The nucleosynthesis at these conditions resembles 
the n-process of Blake & Schramm (1976), but there it was 
discussed to occur because of a decrease of the neutron den- 
sity below 10 18 , while here it happens because of a decrease 
of the temperature and an associated strong reduction of 
the (7,n)-rates. It should be noted that the prominent odd- 
even effect in the isotope distributions in the right panel of 
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Fig. 10. Left: Height of the third abundance peak relative to the second, measured in terms of the ratio (Y (A 3 ) / Y (A2)) , 
as function of the temperature Tg(to) for a power-law decay of the temperature during the second, slow expansion 
phase (blue curve) compared to the dependence of this quantity on Tg(t > to) = Tg in the case of constant asymptotic 
temperature (red line, see also left panel of Fig. [6]). The other wind parameters are again the same as in Figs. 0] [9] 
(s = 145, Y c = 0.42, Tdyn = 5 ms). The right panel shows the distribution of cadmium isotopes at time t = 0.1 s for the 
values of Tf(io) listed in the inset; the neutron number densities at t = 0.1s are given, too (cf. Fig. [8]). The plot should 
be compared with the right panel of Fig. [6j 



Fig. [U for Tg = 1.4 and 0.4 has practically disappeared in 
Fig- EH where the isotope distributions for the same values 
of the transition temperature are much smoother. 

We point out here that the r-processing of heavy nu- 
clei through a quasi-equilibrium of (n,7)-reactions and f3- 
decays at conditions where photodisintegrations are prac- 
tically unimportant was recently also discussed by Wanajo 
(2007), who coined the term "cold r-process"0. We prefer 
to call it "r/3-process" , because this name is conform with 
the denotation of other processes (r-process, rp-process, up- 
process,...) and reflects the essential aspect that character- 
izes this variant of the rapid neutron-capture process. 

3.3. Variations of wind parameters in the exponential phase 

In this section we will discuss the sensitivity of our nu- 
cleosynthesis results to variations of the characteristic out- 
flow conditions like exponential expansion timescale, en- 
tropy, and electron fraction. Table Q] lists corresponding 
parameter values for some of the considered outflows: Td yn 
is the exponential expansion timescal^l, Tg(to) and to the 
temperature and time at which the transition occurs from 
the exponential first expansion phase to the second phase 

3 Different from us, Wanajo (2007) used neutrino-driven wind 
trajectories obtained as solutions of the steady-state wind equa- 
tions, and assumed the transition to a constant freeze-out tem- 
perature Tf at some radius. As argued in Sect. 12. ll such a more 
realistic description of the early outflow dynamics in combina- 
tion with constant conditions at late times shares some basic 
features with the simple two-stage expansion behavior consid- 
ered in our work. 

4 We stress that comparing our expansion timescale with those 
given in other papers requires some caution. Otsuki et al. (2000), 
Sumiyoshi et al. (2000), and Terasawa et al. (2002) defined the 
expansion timescale as e-folding time at T = 0.5 MeV, which is 
compatible with our use. In contrast, Wanajo et al. (2001, 2002) 
defined it as the cooling time from T = 0.5 MeV to T = 0.2 MeV, 
which is about 10% smaller. 
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Fig. 11. Time evolution of the temperature for the different 
expansion histories considered in Fig. 1121 All have an expo- 
nential first expansion phase and second stage of power-law 
decline. Note that the lines coincide for cases with different 
entropies s but the same expansion timescale Td yn - 



with power-law decline of temperature and density (the 
temperature evolution for different considered cases is dis- 
played in Fig. fTTj) . i>; ni = -Rini/ r dyn is the initial out- 
flow velocity at an assumed initial radius i?i n ; = 10 km, 
vq = Wini exp^o/idyn) is the outflow velocity at transition 
time tg , (^4) is the average mass number of nucleosynthesis 
yields, and (Y (A3) / Y (A2)) the ratio of the height of the 
third abundance peak relative to the second as defined in 
Sect. O 

As discussed in Sect. [5J the strength of the r-process is 
determined by the abundance of seed nuclei that is able to 
form, and by the number of remaining free neutrons when 
the r-processing sets in at Tg « 2-3. For the homologous 
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Fig. 12. Time evolution of the free neutron density (left) and final abundance distributions of r-process nuclei (right) for 
outflows with different assumed entropies s and exponential expansion timescales Td yn as listed in the inset. Different 
entropy values are associated with different line colors, a variation of the expansion timescale with dashed lines. In all 
cases we assumed Y c = 0.42 and a transition temperature of Tg(to) = 1 between the exponential first expansion phase and 
the second phase of power-law temperature decline. The corresponding temperature evolution is plotted in Fig. 1111 The 
cyan line denotes the most extreme considered case with an entropy of s = 200 and a dynamical timescale of Td yn = 1 ms. 
It serves for demonstrating the influence of fission cycling (see also Fig. W^i . 



outflows considered there we saw that a larger dynamical 
timescale increases the efficiency of seed formation at the 
expense of free neutrons, and for a very short dynamical 
timescale the expansion can be so rapid that a particles 
and neutrons do not have any time to assemble to very 
heavy elements. In both cases a strong third abundance 
peak cannot develop. A smaller value of the outflow entropy 
also allows neutrons, protons, and a-particles to recombine 
to nuclei more efficiently and therefore lower entropies have 
a similar effect as a slow expansion. 

For a given value of the entropy, significant production 
of high-mass elements therefore on the one hand requires 
that the expansion timescale is sufficiently long, corre- 
sponding to a critical lower limit of the seed abundance 
at the start of the r-processing. On the other hand, the ex- 
pansion must be enough fast because otherwise the seed 
production exceeds Y" s cr too much, and the platinum peak 
cannot be formed because of a disfavorably low neutron- 
to-seed ratio. Our systematic runs for homologous outflows 
showed (see Fig. ^ that with an entropy of s — 145 and an 
electron fraction of Y c = 0.42 only exponential timescales 
Tdyn of less than about 10 ms lead to a neutron number frac- 
tion of Y n J; 0.05 and thus to neutron-to-seed ratios around 
20 at the beginning of the r-process. Only then the nuclear 
flow has a chance to go beyond a mass number of 130 and 
to reach the range of A ~ 196, although the third peak may 
still be significantly underabundant compared to solar val- 
ues (see below). The critical limit for the seed abundance 
turned out to be y s cr - 0.003 (cf. Sect. Of course, this 
number depends on s and Y c , and the homologous expan- 
sion timescale that enables third peak formation is shorter 
for lower entropy values. 

Let us now discuss the case in which the homologous 
outflow with its exponential density and temperature de- 
cline is replaced by a slower power-law temperature decay 
during the late expansion stage, i.e. after the freeze-out of 
charged-particle reactions. This makes strong r-processing 



up to the platinum peak possible for a much wider range of 
dynamical timescales Td yn than in the case of purely homol- 
ogous evolution, practically for all values below some upper 
limit. This can be seen in Fig. [T2l which displays the time- 
evolution of the neutron number density and the final abun- 
dance distribution for different choices of entropies and dy- 
namical timescales. The second stage of power-law temper- 
ature decline is assumed to set in at t — to with a transition 
temperature of X| (to) = I ■ A pronounced platinum peak 
develops for s = 105 if Td yn 2.5 ms, for s = 145 this needs 
Td yn 5 ms, and for s ~ 170-200 it requires Td yn 10 ms. 
An extreme case with s = 200 and Td yn = 1 ms leads to 
fission cycling and demonstrates that even for very fast ho- 
mologous expansion during the first stage the slower evolu- 
tion in the second phase allows all neutrons to be captured 
into heavy nuclei. The dashed lines belong to cases where 
the expansion in the homologous phase is too slow for third 
peak formation. 

Figure demonstrates that for a wide range of 
timescale-entropy combinations in the considered intervals, 
2.5 ms ^ Td yll ^ 10 ms and 100 55 s ^ 200, heavy r-process 
elements up to the third abundance peak and beyond can 
be produced. Inspecting our results we find that a very 
strong platinum peak appears for conditions that roughly 
fulfill the relation s <; 10(rd yn + 10). 

In Fig.[l3]the sensitivity of the r-processing to variations 
of Y e is shown in the case with s = 170 and Td yn = 5 ms. 
In the left panel one can see that a change of Y c from 0.42 
to 0.46 leads to a reduction of n n at the beginning of the 
r-process and a more rapid exhaustion of free neutrons. 
The yields beyond the second peak are correspondingly 
lower, with a growing discrepancy at higher mass numbers 
A (right panel of Fig. [T3J While for Y c = 0.42 the platinum 
peak is slighly overabundant compared to observations (see 
below), increasing Y by about 10% to 0.46 is enough to lead 
to a significant underproduction. For the considered short 
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Fig. 13. Time-evolution of the free neutron number density (left) and final abundance distribution of r-process nuclei 
(right) for Y c — 0.42 (red lines) compared to Y c = 0.46 (blue lines). In both cases we used s = 170, Td yn = 5ms, and a 
transition temperature of Tg(io) = 1 between the exponential first expansion phase and the second phase of power-law 
temperature decline. 



' 10 
: 1 
0.1 
0.01 
0.001 
0.0001 











■ i 






















\<-\ ^^obServed\ ratio 


— > 




















































\ * ^\ 
























\l05 \ 145\ 


\ 170 




i i ill i i i i 


. \ i 



6 8 10 



< 

EN 

< 

g 

< 

v 



10 



0.1 



0.01 



0.001 



<Y(A 3 )/Y(A 2 )> 



E / 
/ 

t 



J Z L 



0.1 



dyn 




Fig. 14. Left: Height of the third abundance peak relative to the second as function of the exponential expansion 
timescale 7d yn , for different values of the entropy (labels at curves). The solid lines correspond to outflows with a power- 
law temperature decay following the exponential expansion, while the dashed line shows for comparison the results for 
s = 145 and constant temperature Tg(t > to) = Tg during the second, slow expansion phase. In all cases the outflow 
was assumed to have an electron fraction of Y e = 0.42 and the transition temperature was taken to be Tg(to) = 1. The 
band between the two bold horizontal lines marks the observational range. Right: Time evolution of the ratio of the 
third abundance peak to the second peak and of the second and third abundance peaks, Y(A 2 ) and Y(A^) (according 
to the definition in Eq. [TTj) . normalized to the total amounts of seed nuclei, l^eed) at any given time t (dashed lines) 
and to the seed abundance at the end of the evolution at t — 2 s (solid lines) , respectively; the different line styles and 
colors are labelled in the inset box. Note that we consider all nuclei with Z > 2 as seeds. The results were computed for 
s = 200, Tdyn — 1 rns, Y c — 0.42, and Tg(to) — 1. The waves visible in the time evolution of the different quantities reflect 
the effects of fission cycling in an environment with a very high free neutron density, which is the case for the extreme 
parameters of this run. 



dynamical timescale this can be compensated by a roughly 
30% higher value of the outflow entropy. 

In Fig. FT?] we provide an overview of the platinum peak 
formation in dependence of the entropy and dynamical 
timescale of the model outflows, using Y c = 0.42 and as- 
suming a transition temperature of Tg(to) = 1 between the 
exponential first cooling phase and the second phase with 
power-law decrease of temperature and density. The left 



panel confirms what we described above: for higher values 
of the entropy the A ~ 196 peak can be assembled for 
an increasingly wider range of expansion timescales. The 
bold horizontal lines in the left panel mark the observa- 
tional band for the abundance ratio (Y (A 3 ) / Y (A 2 )) . The 
third peak relative to the second tends to become over- 
abundant compared to the observations when the expan- 
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sion timescale is low, whereas it remains too weak for long 
expansion timescales. 

The dashed line in the left panel of Fig. [Ml corresponds 
to runs with s = 145 where a constant density and tem- 
perature were adopted during the second phase instead of 
the power-law behavior. The differences between the dashed 
and solid lines for the same entropy are large, in particular 
for longer dynamical timescales. This demonstrates again 
the importance of the late-time behavior of the outflow. 
This importance, however, is significantly reduced when the 
transition temperature Tg(to) is chosen to be near 0.4 (see 
the left panel of Fig. [TO]). 

For entropies near 200 or higher fission cycling can oc- 
cur. The right panel of Fig. [TJ] shows the time evolution 
of the second and third abundance peaks and of the ra- 
tio between both for a case with s = 200, 7d yn = 1 ms, 
and Y c = 0.42. The periodic changes of the abundances 
and of (Y(A 3 )/Y(A2)) reflect the repeated propagation of 
nucleosynthesis wave through the transuranium region via 
fission cycling (the number of nuclear species is increased 
roughly by a factor of three due to fission), similar to what 
was observed by Panov et al. (2003), Goriely et al. (2005), 
and Bcun et al. (2008). The overabundance of the platinum 
peak in this case depends significantly on the nuclear data 
(jPanov et al. 200Tb|) and should be explored separately on 
the basis of new fission data calculations (Panov et al. 
2005). 



4. Discussion and conclusions 

Focussing on the high-entropy neutrino-driven wind envi- 
ronment, we investigated the sensitivity of r-process nu- 
cleosynthesis to the asymptotic low-temperature behavior 
that is assumed to follow an approximately homologous 
early expansion of very rapidly accelerated outflows. On 
the one hand we compared the continuously homologously 
evolving case to a situation where the temperature and 
density reach a lower limit and remain constant later on. 
The latter setting is equivalent to the choice of a freeze- 
out temperature for supersonic wind solutions by Wanajo 
et al. (2002) and Wanajo (2007), and it is similar to the 
use of a constant value for the outer boundary pressure of 
subsonic breeze solutions as adopted by Sumiyoshi et al. 
(2000), Terasawa et al. (2002), and Wanajo et al. (2001). 
Breeze solutions, however, are characterized by the very 
important feature that they imply a causal connection be- 
tween the physical conditions at the neutron star and at the 
outer boundary, i.e. the mass-loss rate and asymptotic tem- 
perature are not independent, whereas supersonic outflows 
with a limiting temperature are not subject to such a con- 
straint. The constant conditions at late times correspond 
to a situation where the homologous outflow is asymptot- 
ically decelerated to very small (zero) velocity. This can 
be considered as a reasonably good representation of the 
situation when a supersonic wind is slowed down in a ter- 
mination shock at colliding with the preceding supernova 
ejecta. On the other hand we investigated a case where the 
velocity was assumed to make a transition from the linear 
radial growth of the homologous phase to an asymptoti- 
cally constant, nonvanishing value. In this case the density 
was assumed to decrease with time like t~ 2 and the tem- 
perature (for constant radiation entropy) like i~ 2 / 3 . This is 
supposed to be a model for the situation when a supersonic 



wind reaches a second stage of slower acceleration at late 
times. 

We have found that for given and constant wind radi- 
ation entropy (we explored s ~ 100-200 fee per nucleon), 
electron-to-baryon ratio (1^ ~ 0.42-0.46), and (sufficiently 
small) wind expansion timescale (rd yn ~ 1-10 ms), a strong 
r-process with production of the third abundance peak de- 
pends not only on the value of the transition tempera- 
ture between the two expansion phases, but also on the 
evolution of temperature and density during the second, 
slower stage. In the case of a constant asymptotic temper- 
ature and density, the formation of the platinum peak is 
enabled when the asymptotic temperature value is moder- 
ately high (T 9 ° - 0.2-0.8). When leaving this range towards 
lower or higher asymptotic temperatures, the possibility 
of third-peak formation strongly decreases. This behavior 
was also seen by Wanajo et al. (2002) and Terasawa et al. 
(2002) for, as they called them, "freeze-out temperatures" 
or "boundary temperatures" above Tg ~ 0.4, but the vari- 
ation with even lower temperatures remained unexplored 
in both works. 

In our second considered case with power-law decline 
of T(t) and pit), the strength of the A ~ 195 abundance 
peak compared to the A ~ 130 peak turned out to be rel- 
atively independent of the transition temperature T f (£o) 
in the interval between about 3 x 10 8 K and 1.4 x 10 9 K. 
For lower transition temperatures the third-peak formation 
decreases steeply and behaves similar to the case with con- 
stant asymptotic temperature and density. This means that 
the appearance and disappearance of a prominent platinum 
peak when the constant asymptotic temperature is lowered 
from Tg ~ 1 to Tg ~ 0.1, is not observed in the case of 
a power-law decline of temperature and density during the 
late wind evolution. Interestingly both considered cases of 
late-time expansion behavior also lead to distinctive dif- 
ferences in the abundances of neighboring isotopes of r- 
process elements. For the power-law time-dependence not 
only more neutron-rich isotopes are formed but the isotopic 
distribution is also smoother. 

Strong r-processing naturally requires a sufficiently fast 
expansion in the homologous phase so that charged-particlc 
reactions freeze out before excessive seed production oc- 
curs and high neutron-to-seed ratios cannot be reached. 
Moreover, a relatively rapid temperature decline during 
the second, slower expansion phase, and at the same time 
a persistence of sufficiently high neutron number densities 
(n n J; 10 24 cm -3 ) are favorable for driving the nuclear flow 
beyond the second abundance peak. For fixed radiation en- 
tropy s as assumed in our models, temperature and density 
are coupled by the relation s cx p/T 3 = const. In this case 
the asymptotic power-law decrease of both quantities en- 
ables third peak formation for a wider range of transition 
temperatures. 

Terasawa et al. (2002) argued that lower asymptotic 
temperatures reduce the charged-particle reactions, which 
leads to less seed production and a higher neutron-to-seed 
ratio, thus causing a better agreement of the r-process 
yields with the solar abundances. For our calculations, 
which are based on the idealized two-stage model of the 
wind expansion, this argument cannot be made. Long be- 
fore the corresponding asymptotic temperatures (between 
Tg = 0.1 and Tg = 1.4) are reached, namely already above 
Tg « 2, charged-particle reactions become inefficient be- 
cause of the impenetrability of the nuclear Coulomb barrier 
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for low-energy thermal protons. Instead, the described nu- 
cleosynthesis results can be understood by the density and 
temperature dependence of neutron-captures and nuclear 
photodisintegration reactions. 

Since the (7,n)-rates decrease steeply with falling tem- 
perature, their role diminishes with a lower transition tem- 
perature. For sufficiently high neutron densities the still 
rapid neutron captures force the r-process path towards the 
neutron-drip line, where the /3-decay rates are large. Instead 
of going through a (7,n)-(n,7) equilibrium, the r-processing 
proceeds now as a quasi-equilibrium of (11,7) reactions and 
/3-decays, which suggests the term "r/3-process" . This sita- 
tion in outflows that expand supersonically and cool quickly 
to a few 10 s K was recently also discussed by Wanajo 
(2007), who named it "cold r-process". The formation of 
the third abundance peak is possible in this situation, pro- 
vided the neutron-to-seed ratio is large enough, because the 
/3-decays are fast and allow for a quick assembling of heavy 
nuclei. In contrast, if the exponential expansion makes a 
transition to a constant temperature that is high, (7,n)- 
reactions drive the r-process path towards the valley of sta- 
bility where the /3-decay rates are small and the r-processing 
therefore slows down. In this case the free neutrons are used 
up in forming the second abundance peak before any signif- 
icant third maximum can build up. On the other hand, the 
power-law temperature decline in the slow expansion phase 
leads to an efficient r/3-process for a wide range of transi- 
tion temperatures Tg(io). Because the strength of the third 
abundance peak then depends only on the neutron density, 
it becomes relatively insensitive to the value of Tg(to) be- 
tween about 0.3 and 1.4. If the transition temperature is 
even lower, the free neutron density during the slow expan- 
sion phase of the outflow is not large enough any more to 
support a strong r-processing. 

In previous investigations (Terasawa et al. 2001, 2002; 
Wanajo et al. 2002) the complex influence of the late-time 
outflow behavior was not explored, but it matters when 
the robustness of abundance yields to variations in the 
neutrino-driven wind models is discussed, especially in sit- 
uations where high entropies, low temperatures, and high 
neutron densities lead to a change of the r-process to an 
r/3-process. 

In summary, we therefore conclude that the detailed 
cooling behavior during the late expansion of supernova 
outflows can have important consequences for r-process nu- 
cleosynthesis. In particular, a slow expansion phase with 
decreasing temperature and density following a rapid, su- 
personic initial outflow expansion, is more favorable for a 
strong r-process than the constant boundary conditions as- 
sumed previously (e.g., Wanajo et al. 2001, 2002; Terasawa 
et al. 2001, 2002). Despite differences in the initial neutron 
number density by more than two orders of magnitude, 
such conditions can lead to a fairly robust production of 
the platinum peak for a wider range of transition tempera- 
tures (at least for the cases of very fast early expansion and 
modest values of entropy and neutron excess as considered 
in our study). This is caused by a kind of self-regulation 
of the quasi-equilibrium of (7,n)-reactions and /3-decays in 
the r/3-process, which leads to roughly the same number of 
neutron captures per seed nucleus, independent of whether 
the process proceeds quickly for high neutron densities (in 
which case the r-process path is closer to the neutron-drip 
line) or slower in the case of low neutron densities. This 
might give a hint for the outflow dynamics of neutrino- 



driven wind ejecta that can enable the uniform production 
of heavy r-process elements suggested by the abundance 
patterns observed in metal-poor stars (see e.g., Cowan & 
Sneden 2006). 

For the explored range of parameter values, the r- 
processing can last 150-200 ms or even longer. This dura- 
tion is the lower time limit for the formation of heavy ele- 
ments from iron-group seeds up to the platinum peak. For 
entropies s between 100 and 200 fee per nucleon, a solar- 
like production of the heaviest elements can occur only 
in the case that the outflow expansion decelerates and a 
transition from the initial homologous phase to a slower 
later stage of expansion occurs. In addition, a short or very 
short exponential timescale (rd yn ~ 10 ms for s = 200 
and Tdyn ~ 1 ms f° r s = 100) is needed during the ho- 
mologous expansion, if moderately neutron-rich conditions 
(Y e — 0.42) are considered. Of course, if the neutron ex- 
cess decreases, the formation of the third abundance peak 
requires a higher entropy. These parameter constraints are 
rather similar to those shown in Fig. 10 of Hoffman et al. 
(1997). This can be understood from the fact that super- 
sonic solutions of the neutrino-driven wind equations as 
those employed by Hoffman et al. (1997) show approxi- 
mately homologous expansion (v oc r) only up to some max- 
imum radius, above which the velocity continues to grow 
significantly less rapidly. Although there is no deceleration 
as it is caused by a wind termination shock (see Arcones et 
al. 2007), the slow-down of the acceleration of supersonic 
winds has still some basic similarity to the transition from 
a rapid first expansion phase to a slower second stage as 
assumed in our simplified outflow models. 
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